Local-field study of phase conjugation in metallic quantum wells with probe fields of both 
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The phase conjugated response from nonmagnetic multi-level metallic quantum wells is analyzed and an 
essentially complete analytical solution is presented and discussed. The description is based on a semi-classical 
local-field theory for degenerate four-wave mixing in mesoscopic interaction volumes of condensed media de- 
veloped by the present authors [T Andersen and O. Keller, Phys. Scripta 58, 132 (1998)]. The analytical 
solution is supplemented by a numerical analysis of the phase conjugated response from a two-level quantum 
well in the case where one level is below the Fermi level and the other level is above. This is the simplest config- 
uration of a quantum well phase conjugator in which the light-matter interaction can be tuned to resonance. The 
phase conjugated response is examined in the case where all the scattering takes place in one plane, and linearly 
polarized light is used in the mixing. In the numerical work we study a two-monolayer thick copper quantum 
well using the infinite barrier model potential. Our results show that the phase conjugated response from such 
a quantum-well system is highly dependent on the spatial dispersion of the matter response. The resonances 
showing up in the numerical results are analytically identified from the expressions for the linear and nonlinear 
response tensors. In addition to the general discussion of the phase conjugated response with varying frequency 
and parallel component of the wavevector, we present the phase conjugated response in the special case where 
the light is in resonance with the interband transition. 



C^ 



X3 

o 



> 

oo 



0\ 

-I— > 



I 

o 
o 



X 



I. INTRODUCTION 

Since the birth of nonlinear opticsEl as a discipline in 
physics nonlinear optical processes have been of great in- 
terest to scientists, for instance tOi-help describe surfaces 
and interfaces of condensed matteroQ Studies of surfaces 
and interfaces of condensed media belong to the regime 
of mesoscopic physics, rWJiere also quantum wells, -wires 
and -dots can be found. QeI Among the many nonlinear op- 
tical phenomena that has been studied in the regime of 
mesoscopicj^hysics are (i) second hamioi^ic generation from 
magneticUliy as well as nonmagnetidlt 
and difference frequency generation,Lj' 
the most Mominent applications, tp^ay is Sisyphus cooling of 
atomsjHirEj (iii) photon drag,EJ"EI 
duced secoodJiaHponic generation 
KeiT effect,BEH"P^ ^i) electronic and vibrational surfaci 



systems, (ii) sum- 
^"' of which one of 



iMX DC-electric-field in- 
|2SI (v) the second-order 



man scattering,EJliJ (vii) two-photon-nbotoemission,! 
(viii) generation of higher harmonics,tj'c2l (ix) the second- 
order Lp£euz.-Mie scattering,cZl and (x) degenerate four-wave 
mixing .c3"EHI 

In the present communication we study phase conjugation 
by degenerate four-wave mixing in a quantum-well structure, 
where both interband and intraband transitions are allowed. 
Phase conjugation is a nonlinear process where the response 
field is counterpropagating to an incoming probe field. The 
usual descriptions (see Refs. p3[-p7| and references therein) 
of degenerate four-wave mixing (DFWM) are based on the 
assumption that the field amplitudes are slowly varying on 
the optical wavelength scale [slowly varying envelope (SVE) 
approximation], and thus also on the electric dipole (ED) 
approximation. We have previously presented the reasons 
(in Refs. M and E9b why these approximations are invalid 
when considering optical interactions with matter of meso- 
scopic size, especially when evanescent components of the 
optical field are present. Four-wave mixing in media with 
two-dimensional translational invariance has so far been stud- 



ied by other authors in the coiitexixif phase conjugation of 
electromagnetic. sjorf ace waves,Eji£]l and of a bulk wave by 
surface waves.EjTj In these investigations macroscopic elec- 
trodynamic approaches were used. In order to go beyond the 
SVE and ED approximations a nonlocal microscopic theoret- 
ical model for optical phase conjugation by DFWM has been 
constructed (see Ref. ^8|) for nonmagnetic media. In addition 
to avoiding the SVE and ED approximations, other usually 
made approximations when considering optical phase conju- 
gation are avoided in our model, namely (i) the paraxial ap- 
proximation, (ii) the assumption of a lossless medium, (iii) 
the assumption of a weak probe field, and (iv) the requirement 
of phase matching between the interacting fields. 

In a two-dimensional translationally invariant system the 
change in energy of an electron due to an electric field can 
either involve a change of momentum along the translation- 
ally invariant plane (intraband transition), a change of energy 
eigenstate perpendicular to the translationally invariant plane 
(interband transition), or both. The change of momentum 
along the translationally invariant plane occurs as an addi- 
tion of the momentum parallel to the plane of the interact- 
ing field component (denoted by q\\ ) to (or subtraction from) 
the momentum of the electron parallel to the surface (de- 
noted K||). It is thus convenient to divide the photon momen- 
tum q into its components parallel and perpendicular to that 
plane, i.e., q = {q\\, q±). Then the vacuum dispersion rela- 
tion q ■ q = q^ = w^/cg provides us with an extra degree of 
freedom, since g[| = \q^^ \ can be larger than oj/cq. Using the 
vacuum dispersion relation we find that q± becomes imagi- 
nary in that case. In the following, propagating field compo- 
nents thus shall refer to the case where q± is a real quantity, 
and evanescent field components to the case where q± is a 
purely imaginary quantity. If we want to get a broad under- 
standing of the phase conjugated response of a probe contain- 
ing both propagating and evanescent field components from a 
quantum- well phase conjugator, two cases are of fundamental 
interest, namely (i) the pure intraband case and (ii) the case 
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where also an interband transition is involved. 

The phase conjugated response from a pure intraband quan- 
tum well we have described in Ref. ^ This analysis re- 
vealed that the phase conjugation reflection coefficient is not 
only highly nonuniform in the q\^ -spectrum, but also that the 
coupling efficiency is several orders of magnitude larger in 
part of the evanescent regime than in the propagating regime. 
Since evanescent waves are strongly decaying in space we 
further concluded that if one wants to see the phase con- 
jugation of evanescent modes, both excitation and observa- 
tion should take place close to the surface of the quantum 
well. Furthermore was discussed the problems of excitation 
of the near-field regime, and the consequences from choos- 
ing a broadbanded (with respect to qji ) two-dimensional point 
source (quantum wire) revealed that parts of the evanescent 
spectrum could be excited, and in Ref. ^ that phase conju- 
gation of evanescent fields can lead to a focus of the phase 
conjugated field substantially below the so-called diffraction 
limit.lI3Q Since this has also been experimentally observed,LJ 
we judge that it is highly relevant also to give an account of 
how evanescent fields are phase conjugated in a system where 
not only one electronic level is present. 

Since including more than one interband transition will 
be necessary for most practical applications, we present in 
this paper the complete solution to the theoretical model of 
Ref. q8| in the case of two-dimensional translational invari- 
ance, although a description based on the_sje]f-field approxi- 
mation according to the Feibelman theoryll3€£l would be suf- 
ficient in order to determine the dominating response. Giving 
a complete solution also allows us to comment on what we 
would lose using the self-field approximation. The solution 
is based on a discretization in the energy levels of the two- 
dimensionally tranlational invariant medium. Contrary to dis- 
cretization schemes performed in real space or Fourier space, 
our discretization does not in itself imply an approximation. 
Thus, once the complete orthonormal set of wavefunctions 
for the phase conjugating medium has been determined, the 
phase conjugated response can in principle be calculated from 
the solution presented in this paper How to find the proper 
set of wavefunctions for a given material system is another 
problem, which for example pfvJbe treated using one of sev- 
eral band-structup|JHethodsJI3'Ea e.g., the Korringa-Kohn- 
Rostoker|4KKR),E3'Ey the linearized augmented plane-wasie 
(LAPW)j!il or the Unear muffin-tin orbital (LMTO) method.H 
These methods are based on an atomic description of the po- 
tential in a certaiiuadius of each atom, adding exchange- and 
correlation termsH and different approximations in the re- 
gions between the atomic boundaries. Using such a method 
one will probably be able to give more accurate numerical re- 
sults for specific materials, but at the cost of the (relatively) 
analytical simplicity. Therefore, we will not elaborate further 
on this point here, but in stead resort to using a simple descrip- 
tion of the matter wavefunctions. Doing so, we will be able 
to present a qualitative discussion based on analytical expres- 
sions. 

Using a two-level quantum-well phase conjugator, it is also 
possible to study resonant four-wave mixing, which until now 
has been studied only without spatial dispersion (at the point 
(g|| , oj) — (0, 0^21) in the gy -a;-plane, UJ21 being the interband 
transition frequency), as described in, e.g., Refs. p4|-p7[ 

Thus, in Sec. O we present the theory in the form of a local- 



field formalism, we choose a scattering geometry, and the so- 
lution is presented as a discretization in the energy eigenstates. 
In Sec. |l| we prepare for a numerical calculation. We start by 
adopting the simple infinite barrier (IB) model to describe the 
quantum well. Furthermore we define the phase conjugation 
reflection coefficient, and the section is concluded with a spe- 
cific choice of a convenient system to investigate. To give 
an impression of the implications of our theoretical model 
we have presented in Sec. IV numerical calculations for a 



two-level quantum- well phase conjugator The calculation is 
supplemented by a discussion of the results, in particular an 
identification of the different resonances appearing when the 
wavenumber along the surface plane as well as the frequency 
varies. In Sec. |v| we widen our discussions, with emphasis 
on (i) the interband resonance, (ii) the validity of the self-field 
approximation, and (iii) the choice of appropriate relaxation 
times. Finally, in Sec. VI, we conclude. 



II. THEORY 

As a forerunner for the analysis of the optical phase conju- 
gation from a two-level quantum well we briefly sketch how 
a local-field calculation allows one to determine the so-called 
degenerate four-wave mixing response of a mesoscopic metal- 
lic film deposited on a dielectric substrate. To create a phase 
conjugated field, which in the plane of the film propagates 
in a direction opposite to that of the probe field, two coun- 
terpropagating pump fields must be present inside the phase 
conjugating medium. Although the theoretical model devel- 
oped in Ref. pq allows us to make almost arbitrary choices 
of the interacting optical fields, we will in the present work 
assume for simplicity that the pump fields (i) propagate par- 
allel to the plane of the film, and (ii) have constant amplitude 
across the film. The scattering geometry is shown in Fig. IT] 
together with the chosen coordinate system. We will further 
limit our study to the case where (iii) scattering takes place in 
the x-z-plane, and (iv) the interacting fields are linearly polar- 
ized, either in (p-polarized) or perpendicular to (s-polarized) 
the scattering plane. Since it is necessary in a study of non- 
linear optical phenomena in mesoscopic interaction volumes 
to abandon macroscopic electrodynamics, the starting point is 
the microscopic Maxwell-Lorentz equations. The phase con- 
jugated field from a quantum well exhibiting free-electron- 
like dynamics in the plane of the well {x-y-plane) can thep be 
described using the single-coordinate (z) loop equatiorO 
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FIG. 1. The system we consider here consists of a three layer struc- 
ture, namely (i) vacuum, extending from — cxj to —d, (ii) quantum 
well, extending from ~d to 0, and (iii) substrate (crosshatched), ex- 
tending from to +CX3. The three incoming electromagnetic fields 
consists of two pump fields (labeled "1" and "2") and a probe field. 
Also shown is the Cartesian coordinate system used in our calcula- 
tions. 
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Epc{z;q\\,uj) == i?pc(z;<7||,cj) - ifj-auj / / G(z, z"; gjj,u;) 

•ct(z",z';(7||,lj) • Epciz';qi\,u!)dz"dz', (1) 

where w is the common angular frequency of the participat- 
ing fields, and qji is the component of the probe field in the 

film plane. It is the so-called background field, Ep(~,{z; gji , uj), 
which makes the loop problem different for the various non- 
linear (and linear) problems. It is here given by 



E 



pc(^;'zih^) 



?(3)/ 



(3), 



i^QUJ I G{z,z';qi\,uj) ■ JlJ{z';qi\,uj)dz', 

(2) 



where J_J (z'; (7J1 , uj) is the current density driving the nonlin- 
ear process. The pseudo- vacuum propagator 0(2, z"; q\\,oj) is 
given by 



G{z,z';q\i,uj) 



^iq±\z-z 



2iq_ 

+Q{z' - z)er (8) eV] + 



g-igj^(z+z') 



2iq± 



J + 9(z - z')e, ®ei 
[r^ey ® Cy ~\- r^Cr ® et 



aij{z,z';q\\,uj) 



2i 1 



2 A^ 

nm 



UJ 



-—Siz 

q^ 



(3) 



where the first term describes the direct propagation of the 
electromagnetic field from a source plane at z' to the observa- 
tion plane at z, the second term accounts for the reflection 
at the quantum-well/substrate interface, and the third term 
characterizes the field generated at the observation plane by 
the current density prevailing in the same plane (thus named 
the self-field term). Above, e^ = q^^{q±,Q,—q^^), and 
eV — q~^{—q±,0,—q\\), taking q\\ = q\\ex. The quantities 
r'^ and r^ are the amplitude reflection coefficients at the vac- 
uum/substrate interface in the absence of the quantum well 
for s- and p-polarized fields, respectively. Both of these are in 
general functions of q\^. Moreover, the vectors Cx, Cy, and e^ 
are unit vectors along the principal axes in the Cartesian x-y- 
z-coordinate system, 0(- • •) is the Heaviside unit step func- 
tion, and S{- ■ •) is the Dirac delta function. 

The Jj'th tensor element of the linear response tensor 
ct(z", z'; qji , Lo), appearing in Eq. (|l]), is given byE^Ea 



fni^W +q\\)- fmiK\\) 



Wnm(K|| + q\\,K\\) Wnm(K|| + q\\,K\\) " W 



jz,nm{z;2Kil + q|| ) jj,„„ (z'; 2k|| +qii)d K||,(4) 



provided the set of wavefunctions is complete. In Eq. (j^ we have introduced the transition current density in the mixed Fourier 
space, namely 



Jnm \Z ^ ' 



2im„ 



^Q^^iJ*^{z)M^) + e. ( C.(^)%^ - V„(z)'^'^"^"^ 



dz 



dz 



(5) 



In relation to Eq. (^, Q|| is equal to 2k|| + qp where K|| is 
the wave vector of the given electron in the plane of the well. 
The transition current density also occurs in the nonlinear re- 
sponse tensor (see Appendix pi) and in this context various 
combinations of g||, fcy, and K|| appear in Q\\. The quantities 
■0a, a € {n,m}, are the one-dimensional electronic energy 
eigenstates of the quantum well belonging to the z-direction, 
and they satisfy the field-unperturbed Schrodinger equation 
Ti-Qipa — Sai^a- The quantity fa{K\\) denote the Fermi- 
Dirac distribution for the eigenstate 'i'a(r) = ^a{z) exp(iK|| ■ 
r)/{2TT), where also the solution to the Schrodinger equation 
along the quantum well is taken into account. It is given by 
/a(K||) = [l+exp{(£a+ft2K|/(2me)-M)/(A:BT)}]-i, where 
/cb is the Boltzmann constant, /i is the chemical potential of 
the electron system, and T the absolute temperature. For the 
various Cartesian components of the transition current density, 
we use the notation ji^nm,{z\ K||), i G {x, y, z}. The complex 
cyclic transition frequency is defined by 



i{Q\\.a,Q\\,b) = 



+ 



2™. 



{\Qla\^-\Q 



\.b\ 



(6) 



where e„ and £„ are the eigenenergies of the quantum well 
states belonging to the z-direction, and Q\\^a and Qn j, can be 
any of the relevant combinations of gj| , fc|| , and K|| . The quan- 
tity Tnm is the relaxation time. 



The nonlinear current density, J_J (z'; qi^\ , uj), is related to 
the pump and probe fields by a constitutive relation of the 
form 



J^-liz;q'\\,u;) 



1 



i27:y 



^(z,z'; (7||,A:||,w) 



:E{-ki\ , Lu)E{k\\ , lu)E* (z'; -g]| , uj)dz' + i.t.. 



(7) 



where 



^(z, z'; (7||, fc||,a-') = / ^{z, z' , z" , z"';q\\,k^^,uj)dz"'dz'' 



(8) 



is the relevant nonlinear response tensor when the pump fields 
are essentially constant (slowly varying) across the quantum 

well, i.e., ii;(z"';-fcj|,w) = £{-^,0;) and E{z";kn,uj) = 

E{k\\ , (jj) in Eq. (^. Within the framework of a single- 
electron random-phase-approximation approach an explicit 
expression for S(z, z', z", z'"; q^\ , fcy , uj) has been established 
in Ref. E^ The term "i.t." denotes the so-called "interchanged 
term", which takes into account the symmetry of the pump 
fields. It is obtained from the first term by interchanging the 
two pump fields (the pump field wavevector fc|| is replaced 
by ~^||)- The explicit expression for the simplified non- 
linear conductivity tensor, S(z, z'; (7J|, fc||, w), can be found 
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in Appendix Q while the more general conductivity tensor 
E{z, z', z", z'"; Oil ,k^^,uj) has been given in Ref. 58, We have, 



however, in Appendix W only listed one of the seven parts, 
namely part G, of the nonlinear conductivity tensor that ap- 
pears in Ref. pq, since when interband transitions are strong, 
it is dominating the response by several orders of magnitude 
compared to the other six (A-F). 

As a consequence of the above-mentioned choice (but inde- 
pendent of the direction in which the pump fields propagate) 
the number of terms in the nonvanishing elements of the non- 
linear response tensor is further reduced, since the orthonor- 
mality of the z-dependent parts of the wavefunction gives 



'4>niz)lpm{z)dz = S„ 



(9) 



where Snm is the Kronecker delta. Also, by integration of the 
microscopic transition current density given by Eq. (g) over z 
one finds 



Jnni{z;Q\\)dz 

where 



eh 
2inie 



^W||^n77i ~^ Pz.nm^z 



(10) 



i^:jz)'-^-Mz)'-^)dz (11) 



dz 



dz 



is proportional to the z-componeoLpf the electric dipole mo- 
ment related to the nm-transition.E3 

The conductivity tensor S(z, z', z", z'"; (tJi, /c[|,cj) has in 
general 8 1 nonzero tensor elements (3x3x3x3) and consists 
of seven different parts (A-G) after the seven different phys- 
ical processes contributing to the response (see Ref. |8] for 
details). When scattering takes place in the x-z-plane with 
linearly polarized light the general treatment can be split into 
eight separate parts related to the possible combinations of po- 
larization of the three different incident fields. In this scatter- 
ing geometry gj| and fc|| lie along the x-axis, giving a mirror 
plane at y — 0. Consequently, only tensor elements in the 
nonlinar response tensor with a Cartesian index even num- 
bered in y contributes, and the 81 tensor elements generally 
appearing are reduced to 41. The separation of the tensor el- 
ements into the eight sets of elements contributing in these 
configurations follows in a straight forward manner from the 

definition of the sum-product operator ":" between the non- 
linear current density and the interacting electric fields, i.e., 

[E:EEE*\, = J2jkh ^ijkhEhEkE*. The added restriction of 
letting the pump fields travel along the a;-axis then reduces the 
number of contributing matrix elements from 41 to 18, since 
when traveling along the x-axis, the pump fields are polarized 
in either the y-direction or the z-direction. The resulting sets 
of tensor elements we have presented in Table |. 

To solve Eq. ([l|), we can establish a so-called coupled an- 
tenna loop. First, we notice that each matrix element of 
the linear conductivity tensor [Eq. (Q)] with the insertion of 
Eq. (^ can be written as a product of a z-independent term 
and two terms depending on z and z', respectively. Element 
ij then takes the form 

(Ty(z,z';gj|,w) =^Q]i^{q\\,Oj)jz,nm{z)jj,mn{z'), (12) 



TABLE I. Contributing tensor elements of the nonlinear conduc- 
tivity tensor when the pump fields are propagating in the x-direction 
and all fields are polarized in (p) or perpendicular to (s) the 
x-z-plane. The left column shows the polarization combination of 
the incoming fields (pump 1, pump 2, probe), the center column 
shows the polarization of the phase conjugated field, and the right 
column shows the tensor elements contributing to the nonlinear in- 
teraction. 
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where ]nm{z) = ]nm{z]ex + By)- The various Q-quantities 
can readily be identified from Eq. (^, and the integrals can be 
solved using the method described in Appendix 0. Inserting 
Eq. ([l|) into Eq. ([I|), we get 

^PC(^) = ^pc(^) + ^^mn(2) • f„„, (13) 



omitting the reference to q\\ and lo for brevity. In Eq. (ol) we 

have introduced the 3x3 tensor F„m(z) with the nonzero 
elements 

i^„-;^(z) - -i^^oc^ ^ Qr4y'G.,(z,z")j.,„™(z")dz" 

qj- 



i\\ 



. p2=^ (z) 



(14) 



Kmiz) == -iA*OW Yl 2'"" / Gxiiz,z")j,^nrniz")dz" 



—F^^ (z) 
9|| 



and the elements of the vector T^n are written 



(15) 



F^liz) = -ifiou^Qll, / Gyy{z, z")jy^^„Xz")dz" , (16) 



^i,,nn= I ji,mn{z')Epcs{z')dz' , iG{x,y,z}. (17) 



To determine the phase conjugated field the quantity Tmn 
must be calculated. This is done by multiplication of each 
element Epc.i{z'), i G {x, y, z} of the phase conjugated field 
in Eq. (|l|) by the relevant ji^mniz) followed by an integra- 
tion over the z-coordinate. Hence, when the phase conjugated 
light is s-polarized, Eq. ( 13 ) is transformed into the following 
set of linear algebraic equations: 



-L j/,mn / y ^yy,mn^ yM — "a,mnj 



(18) 



i.e., n? equations with just as many unknowns. In the case of 
p-polarized light, we obtain 
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^ x,mn / ^ \^ XX ,mn'^ X ,vl ' -'^xz.mn^ z,vl ) "a;,mn7 \^^) 

vl 

J- z,mn ~ / ^ \^zx,mn'^ x.vl + ^zz.mn'^ z.vl) ~ "z,mnj l^u) 



vl 



which are 2n^ equations with just as many unknowns. In 
Eqs. (p^-(po[) above, the elements of the vectorial quantity 
(imn are given by 



^J,mn = I ji,rnn{z)Ep(.j^{z)dz, ie{x,y,z}, (21) 



and the 3x3 tensorial quantity K 
nonzero elements 



vl / 
mn \ 



, oj) has the five 



K. 



vl 



ji,mn{z)Fil{z)dz, 



(22) 



where the indices "i and "j" can take the values of ij e 
{xx,xz,yy, zx, zz}. By means of the procedure sketched 
above, we have been able to transform the integral-equation 
problem for the phase conjugated field, Epc{z), [Eq. (|])] to 
a matrix problem for the Fmn-vectors. This discretization in 
the energy levels is exact, and once the linear algebraic set of 
equations for the Fmn-vectors, truncated so as to keep only 
the subspace of relevant energy levels, has been solved (nu- 
merically) the phase conjugated field can be obtained from 
Eq. (|3]). Integral equations of the type given in Eq. ([l|) is 
often solved (numerically) by discretization in the real space 
coordinate. By such a procedure one has to worry about how 
small discretization lengths one may dare to take from a phys- 
ical point of view. The discretization in energy levels used 
here does not suffer from this uncertainty. 



III. NUMERICAL FRAMEWORK 

Our description of the phase conjugated field has until now 
been independent of the actual wavefunctions in the quantum 
well, and thus also independent of the form the potential takes 
across the active medium. However, if we want to perform a 
numerical calculation of the phase conjugated field we have to 
choose a definite potential across the quantum well, giving us 
a set of wavefunctions to work with. Below we use the infinite 
barrier (IB) model potential for the numerical study, since this 
model is sufficient for a qualitative study. 

As shown in Fig. 0, in this model the one-dimensional 
potential V{z) is zero inside the quantum well (in the in- 
terval —d < z < 0) and infinite everywhere else. The 
stationary state wavefunctions inside the quantum well are 
given by tjjn{z) = y^2/dsm{nnz/d) and outside the quan- 
tum well, ipn{z) = 0. The associated eigenenergies are 
£„ = {mrhy/{2m,,d^). Within the IB model, Eq. ^ gives 



Pz 



4nm[l - (-1)"+™] 
(n^ — im?)d 



(23) 



for n ^ m, and Pz.nm = for n — m. For a metallic quan- 
tum well one may even at room temperature approximate the 
Fermi-Dirac distribution functions by their value at zero tem- 
perature, i.e.. 
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FIG. 2. Infinite barrier (IB) model potential (thick solid line) for a 
quantum well with boundaries at 2 = —d and 2 = 0. In the present 
case, only one energy level below the Fermi energy (here called |1), 
with energy ei) and one energy level above the Fermi energy (called 
1 2), with energy £2) contributes to the solution. The remaining in- 
finite set of energies appearing in the IB model we assume are so 
far away from |1) and |2) that they do not contribute to the solution. 
The dotted curves indicate the shape of the wave function for each 
of the two energies. To the right is shown the possible transitions, 
identified with their respective transition frequency and relaxation 
time. 
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where £p is the Fermi energy of the system. In the low- 
temperature limit it is possible to find analytical solutions to 
the integrals over K|| appearing in Eq. ( Al ). The explicit cal- 
culations are tedious but trivial to carry out, and since the final 
expressions are rather long we do not present them here. For 
the interested reader some steps in the calculations are repro- 
duced in Appendix g[ 

The Fermi enejEgy is calculated from the global charge neu- 
trality condition,E3 aad for a quantum well described by the 
IB model, it becomesEa 



irh^ 



NpTUf, 



ZN+d- 



Nf{Nf + 1){2Nf + ly 



2d2 



(25) 



where N+ is the number of positive ions per unit volume, Z 
is the valence of these ions, and Nf is the quantum index of 
the highest occupied level. From Eq. (|5]), the number of oc- 
cupied levels can be calculated if the thickness is known, and 
vice versa. The minimal thickness for the quantum well to 
have n levels below the Fermi level can be determined from 
the relation £f = £n, and the maximal thickness from the 
condition £f ~ Sn+i- Thus for n bound states below the 
Fermi energy we find the minimal and maximal thicknesses 






2ZN, 



(n+l)(2n+l) 
6 



(26) 



i.e., a result that depends on the number of levels below the 
Fermi energy and the number of conduction electrons in the 
film. 
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To estimate the amount of phase conjugated light, we use 
the phase conjugation (energy) reflection coefficient defined 
as 






Iprohc{~d;q\\,uj)' 



(27) 



in which I^^\ I^^\ /probe, and /pc are the intensities of 
the two pump beams, the probe and the phase conjugated 
field, respectively. Each of the intensities are given by / = 
^eaCoE ■ E*{2TTy^. Thefactorof(27r)"'' originates from the 
manner in which we have introduced the Fourier amplitudes 
of the fields. 

For the remaining part of this work we choose a copper 
quantum well with N+ = 8.47 x lO^^m"^ and Z = 1 (data 
taken from Ref. pTl). The Cu quantum well is assumed to be 
deposited on a glass substrate for which we use a refractive in- 
dex 71 of 1 .5 1 . With this substrate, the linear vaccum/substrate 
amplitude reflection coefficients can be obtained by use of the 
classical Fresnel formulae r" = [qj_ — [n^q^ — q},Y/'^]/[q±^ + 

(n2q2 _ gpl/2] ^„j ^p ^ [jj2^^ _ (^^2^2 _ ^2)l/2]/[„2^_^ ^ 

{v?q^ — (7y )^/^]. Having the pump wavevectors parallel to the 
X-axis then gives a pump wavenumberof fcii ~ nq = 1.51q. 



IV. NUMERICAL RESULTS FOR A TWO-LEVEL 
QUANTUM WELL 

To calculate the phase conjugated response from a quan- 
tum well with an arbitrary number of bound eigenstates one 
would have to superimpose interband and intraband contribu- 
tions. Thus in a study of the complete response where local- 
field effects are neglected one basically would have to add 
the contributions from the various pairs of levels located in 
different subbands or in the same band. Seen in this light, 
thorough treatments of the single-level case, where only in- 
traband transitions are allowed, and the two-level case, where 
transitions between two eigenstates located in different bands 
occur, would form a good qualitative starting point for anal- 
yses of multi-level quantum-well systems. The single-level 
case we have studied beforejlj and the following treatment 
will thus be directed towards a description of the phase con- 
jugated response from a two-level quantum well. Thus, we 
choose the simplest possible configuration in which interband 
transitions can occur, i.e., a quantum well with only one bound 
state below the Fermi energy. Above the Fermi energy we 
also assume that only one bound state can be reached, and 
thus the wavefunctions are ipi{z) = ^y2/dsm{nz/d) and 
■02(2) = ^/2/dsin{2^Tz/d). The associated energies then 
become ei = (7r/i)^/(2med^), and £2 = (27rft)^/(2med^), 
respectively. The quantum well with the various relevant en- 
ergies and wavefunctions, as well as the electronic excitations 
are shown in schematic form in Fig. g. In the present two-level 
case Eq. (EJ) becomes 



16 

Pz,n7n = -T-:Sgn{n ~m), 
6a 



(28) 



where {n,m) G {(1, 2), (2, 1)}. If just the ground state 
should have an energy less than the Fermi energy, we see from 
Eq. (pq) that the film thickness must be less than dmax — 



^3Tr/{2ZN+). The minimal thickness is in the IB model 
zero, but in reality the smallest thickness is a single mono- 
layer Using Eq. (|6|) the maximal thickness for a two-level 
Cu quantum well then becomes dmax ~ 3.82A, which is more 
than two monolayers and less than three. Thus we have two 
obvious choices for the thickness of the quantum well, namely 
a single monolayer or two monolayers. We choose two mono- 
layers, since by this choice the two energies ei and £2 are 
closest to each other, and thus the energy needed for a reso- 
nant transition to occur is lowest. Two monolayers of copper 
roughly corresponds to a thickness of d = 3.6 A (bulk value). 
With this choice, the energy difference between the two states 
is £2 — £1 = 8.70 eV, and the corresponding resonance in the 
optical spectrum is found at the wavelength A — 142.4 nm. 



A. Phase conjugation reflection coefficient 

Among the eight possible ways of using linearly polarized 
light in our chosen scattering configuration, two combinations 
give an s-polarized response when using an s-polarized probe 
field, the pump fields being either s-polarized or p-polarized, 
but with the same polarization for both pump fields. When the 
pump fields are s-polarized, the nonlinear conductivity tensor 
element that contributes to the response is 'E.yyyy. Altogether 
the phase conjugated response in this purely s-polarized case 
(called "sss") is negligible, since it is tens of orders of mag- 
nitude less than those of the other combinations. If, on the 
other hand, the pump fields are p-polarized (pps), '^yyzz is flie 
element of the nonlinear conductivity tensor that contributes. 
Plotted as isophotes (contours of equal intensity) in the nor- 
malized tj-g|| -plane (uj normalized to the interband transition 
frequency u!2i and q\\ normalized to the vacuum wavenumber 
uj/cq), the result is shown in Fig. ^ 

Two other combinations of polarization give p-polarized re- 
sponse using a p-polarized probe field. As above, the pump 
fields have to be of the same polarization, and can either be s- 
or p-polarized. With s-polarized pump fields (ssp), four ten- 
sor elements of the nonlinear conductivity tensor contribute 
to the phase conjugated response (see Table |b. The phase 
conjugated response is shown in the normalized tj-g|| -plane in 
Fig. H In the other case, another four tensor elements of the 
nonlinear conductivity tensor contribute to the phase conju- 
gated response when the pump fields are p-polarized (see Ta- 
ble ||). We have in Fig. g shown the phase conjugated response 
for this configuration (ppp) in the normalized tj-gn -plane. 

In the remaining four cases, the response has a different po- 
larization than the probe field. This is obtainable by the use 
of differently polarized pump fields. In order to achieve an s- 
polarized response from a p-polarized probe field one makes 
use of two differently polarized pump fields, and four tensor 
elements of the DFWM response tensor contribute to the so- 
lution, cf. Table |. Similarly, two differently polarized pump 
fields are needed in order to produce a p-polarized response 
from an s-polarized source. For this process, another four ten- 
sor elements of the nonlinear conductivity tensor contributes 
according to Table |. Since the resonance structure of these 
last four cases are similar, it is sufficient here to discuss the 
result obtained for just one of those cases. Thus, in Fig. the 
result is shown for the case where pump field 1 is s-polarized 
and pump field 2 and the probe are p-polarized (spp). 
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FIG. 3. The phase conjugation reflection coefficient from a 
two-level metallic quantum well is plotted in the case where 
s-polarized probe field gives s-polarized phase conjugated response, 
and where the pump fields are p-polarized (pps). The response 
is plotted as isophotes (contours of equal intensity) [m'^AV'^] on a 
logarithmic scale as a function of (i) the frequency uj normalized 
to the transition frequency lu2i, and (ii) the parallel component of 
the wavevector, normalized to the vacuum wavenumber. The dif- 
ference between two neighbouring contours is one order of magni- 
tude. To indicate the absolute amplitude, the isophote of value 10^'^° 
m^AV^ has been plotted using a long-dashed curve and the isophote 
with magnitude 10^'*" m^AV^ with a short-dashed curve. On the 
Qll/g-scale, the response has been plotted on a linear scale in the 
range < q\\ /g < 0.1 and on a logarithmic scale above qy /q = 0.1. 




FIG. 4. The phase conjugation reflection coefficient from a 
two-level metallic quantum well is plotted in one of the cases where 
p-polarized probe field gives p-polarized phase conjugated response. 
In this case the pump fields are s-polarized (thus named "ssp"). The 
response is plotted as isophotes [m^AV^] on a logarithmic scale as a 
function of (i) the frequency u) normalized to the transition frequency 
uli\2, and (ii) the parallel component of the wavevector, normalized 
to the vacuum wavenumber. The difference between two neighbour- 
ing isophotes is one order of magnitude. Again, the two isophotes of 
magnitude 10"^" m'^AV^ and 10"^° m'^AV^ has been plotted with 
long- and short-dashed curves, respectively. As before, below 0.1, 
§11 /q has been plotted on a linear scale while above it is logarithmic. 
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FIG. 5. The phase conjugation reflection coefficient from a 
two-level metallic quantum well is plotted in the other case where 
p-polarized probe field gives p-polarized phase conjugated response, 
this time with p-polarized pump fields (ppp). As in Figs. H and n, the 
response is plotted as isophotes [m*AV^] on a logarithmic scale as a 
function of (i) the frequency uj normalized to the transition frequency 
0J12, and (ii) the parallel component of the wavevector, normalized to 
the vacuum wavenumber. Again, the difference between two neigh- 
bouring contours is one order of magnitude, and as before, the long- 
and short-dashed curves represents magnitudes of 10^ ' m AV and 
10~^° m*AV^, respectively. In the big picture, q\\/q is plotted on a 
logarithmic scale, while in the strip it is plotted on a linear scale. 



FIG. 6. The phase conjugation reflection coefficient from a 
two-level metallic quantum well is plotted in one of the cases where 
p-polarized probe field gives s-polarized phase conjugated response. 
In this case, pump field 1 is s-polarized while pump field 2 is 
p-polarized (spp). The response is plotted as isophotes [m^AV^] on 
a logarithmic scale as a function of (i) the frequency uj normalized to 
the transition frequency a;i2, and (ii) the parallel component of the 
wavevector, normalized to the vacuum wavenumber. The difference 
between two neighbouring contours is one order of magnitude. The 
absolute amplitude of the isophote of value 10~^" m^fW^ has been 
plotted using a long-dashed curve and the isophote with magnitude 
10^^" m'^AV'^ with a short-dashed curve. The strip below is plotted 
in a linear scale in gy /q while the rest is on a logarithmic scale. 
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The IB model only offers a crude description of the elec- 
tronic properties of a quantum well, since, for example, the 
electron density profile at the ion/vacuum edge is poorly ac- 
counted for. This gives too sharp a profile and underestimates 
the spill-out of the wavefunction. Altogether one should be 
careful to put too much reality into the IB model when treat- 
ing local-field variations (related to, say, ^n or q±) on the 
atomic length scale. Furthermore, neclecting the Bloch char- 
acter of the wavefunctions accounting for the dynamics in the 
plane of the well is doubtful in investigations of the local field 
among the atoms of the quantum well. The crucial quantity 
in the above-mentioned context is the Fermi wave number 
kp = (2me£^i;')^/^//i, and in relation to Figs. ||-|[ only re- 
sults for q\\/q ratios less than approximately 



kp 
Q 



ZN+d ^ 1 



27r 



(29) 



appears reliable. Thus we have cut off our results at the line 
q\\/q = kp/q in the ti>/aj2i-'7||/'?-plane in Figs. ||-^. 

In many theoretical studies of the properties of phase con- 
jugated. &^lds it is assumed that the phase conjugator is 
ideal.E3~0 By this is meant that the phase conjugation re- 
flection coefficient is independent of the angle of incidense of 
the (propagating) probe field (and maybe also of the state of 
polarization). As we concluded for the single-level quantum 
well,c£l and as we can now see for the two-level quantum well 
in Figs. H-ra this assumption is not such a good approximation, 
at least not for a metallic quantum well system. 



B. Resonant structure of the DFWM reflection coefficient 



Looking at Figs. ^|-^ a number of resonances occur They 
can all be accounted for from the analytic solution to Eq. (^Ij) 
by looking at the denominators appearing in the analytic de- 
composition of the products, as given by Eqs. (B2) and ( [B3| ) 
in Appendix pi These resonances are shown on the scale 
of Figs. g-|] in Fig. 0. In the analytic solution of the inte- 
grals over K|| shown in Appendix 0, the solution to the terms 
with three multiplied denominators is reduced in Eq. (B3) 
to the problem of finding a basic solution to the integrals 
over K|| for each of these denominators multiplied by a K||- 
independent factor. The resulting integrals do not contain 
sharp resonances, but the factors in front of them do, when 
aihj — biGj — 0, for i,j e {1,2,3} and i ^ j. In or- 
der to make an analytical treatment of the resonances ap- 
pearing in the nonlinear conductivity tensor we in the fol- 
lowing define a term of the nonlinear conductivity tensor as a 
product of three denominators in Eq. ( Al), and number them 
1,2,..., 12. However, not all terms gives contributions to the 
result in a two-level quantum well. The terms that does not 
give any contributions are the terms with a 2cj-contribution in 
the denominator, i.e., terms 1-2 and 1 1-12. When the denom- 
inators of the rest of the terms (3-10) are put into the form of 
Eq. ( B 1 ), a total of four different as and nine different 6's ap- 
pear. They are listed in Appendix y. Since we are looking for 
the location of the resonances in the system it is reasonable 
in the following analysis to let the respective relaxation times 
Tnm in Eqs. (C5)-(C13) be infinite. 




10 



lOi 



10° 

UJ/UJ21 

FIG. 7. Resonances of the nonlinear conductivity tensor are 
shown as a function of (i) the optical frequency normalized to the 
transition frequency (u>/u!2i) and (ii) the parallel component of the 
wavevector normalized to the vacuum wavenumber (gp /q). This fig- 
ure shows only the pure resonances. The broadening due to the re- 
laxation times is neglected by setting them all to infinity. 



In terms of the a's and &'s listed in Appendix 0, we observe 
that the third term of 'E.ijkh has resonances at (i) ai6f„„ — 
bl^a2 = 0, (ii) a^bi^^-bl^a, = 0, and (iii) a^bl-bl^a, = 0. 
After insertion of the relevant a's and 6's, substitution of fc|| in 
favor of nuj/co (since fc|| — nuj/co in our treatment), and a 
normalization of (7|| to q, i.e., gn = {q\^/q)uj/co, we may solve 
the resulting second order equations with respect to lj as a 
function of gn /q. Then resonance condition (i) gives 
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niecl n-q\\/q / 



± 



me4 n-q\\/q 



2mccl 



hnq\i /qn + q\i/q \ \ hnq\^ /qn + q\\ /q J h?{n + qii / q) 
resonance condition (ii) becomes 



^m ^n 



q\\/q 



(30) 



f^nq\\/q \ \fmq\i/(, 



2mecl 



Q 



^rn ^n 



n + q\\/q q\\/q 



(31) 



and condition (iii) is 



'e'-O 



± 



hnq\\lq y \hnq\\/q) h?q\\/q 



2mp 



In some of the above equations, some of the solutions can 
be ruled out immediately, since, for example, in Eq. (|T]) the 
minus in front of the square root gives only rise to negative 
values of ui in the "interchanged term" (when A:|| — —nuj/cQ). 
In the fourth term of 'E.ijkh we observe that in addition to 
a resonance of type (ii), resonances appear at (iv) a26^„ — 
bl^vO-2 = and (v) a^b^^ - Ki'^2 = 0. Again inserting the re- 
spective a's and 6's from Appendix 0, substituting ncj/co for 
/c|| , and normalizing gy to the vacuum wavenumber, resonance 
condition (iv) becomes 



1 
2h 



y€y ^m}-! 



(33) 



and resonance condition (v) is equivalent to Eq. ([3l|), taking 
into account the interchanged term. In our configuration, the 
choice of a two-level quantum well puts some restrictions on 
the values of the quantum numbers n, m, v, and / in order to 
get a nonzero result. 



Comparing Eqs. ( |A1[ ) and ( |10| ) we ob- 
serve that if pump field one (indexed k) is s-polarized then 
I = V, while / 7^ w if it is p-polarized. Similarly, if the other 
pump field (indexed h) is s-polarized we get m — I, while we 




(32) 



get m 7^ Zif it is p-polarized. These conditions are summa- 
rized in Table ||, and the contributions from Eqs. (p0|)-(p3|) to 
the resonances in Fig. are shown in Table [Il| for the valid 
combinations of quantum numbers. 

TABLE II. Restrictions on the valid combinations of quantum 
numbers for a two-level quantum well in the nonlinear conductiv- 
ity tensor for the three combinations of polarized light of the pump 
fields treated in this communication. Pump field 1 is indexed k, and 
pump field 2 ft in Eq. (41 1. 



k h 



s s 
s p 

p p 



^ terms 3^ 



^ terms 5-8 



^ terms 9-10 



I — V = m 
I — V Am ^ I 
m^ I Al ^v 



V — n Am ^ I 

V = n Am ^ I 



V = n — I 

V = n Al =/= V 



V ^ n Am, ^ I V ^ n Al ^ V 



The resonances conditions in the fifth term of Eij^h are (vi) 
aibtrn + bLo'2 = 0, (vii) a4fof„„ - 6^,„a2 = 0, and (viii) 



'^4^im + ^vm'^1 = 0- ^y insertion of the respective a's and 6's 
from Appendix 0, substitution of fc|| by uuj/cq, and normal- 
ization of q|| to the vacuum wavenumber, resonance condition 
(vi) becomes 






± 






27TT.eCgn 



r9{n-q\\/q)q\\/q 



q\\/q 



f-v -e; 



(34) 



condition (vii) gives 



hnq\\lq 



± 



\hnq\\/q 






q\\/q 



■q\\/q. 



(35) 



and case (viii) becomes 



OJ 



ruecl 



± 



( rUecl 



hnq\\/q y \hnq\\/qj hPq^\/q 



2meCo 



n-q\\/q. 



(36) 



The sixth term of ^yfeh has aresonanceof the type (vi), and further resonances at (ix) a2fo*,,„— 6^;a2 — Oand(x) aib'[i — h^^a2 = 
0. Insertion of the different a's and 6's, k\\ — nu/cQ, and nomalizing q\\ to the the vacuum wavenumber gives (ix) resonances at 



hnq\\/q 



± 



f'e^{) 



hnq\\/qj hP'nq^Jq 



^i\ 



(37) 



and (x) resonances equivalent to those given in Eq. 



In the seventh term of :iijkh there is a resonances of the type of case 



(ix), and furthermore at (xi) a^bt„i — b^iO-2 — and (xii) 03 6^^ — ^ni^2 — 0. As in the previous cases we insert the different 
a's and 6's found in Appendix y, replace /c|| with uuj/cq, and nomalize q\\ to the the vacuum wavenumber. Then case (xi) gives 
resonances at 
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'■e'^O 



q\\/ci 



hnq\i/qn-q\i/q 



± 



(■eCg 



\hnq\\/qn- 



-q\\/q 



q\\/q 



+ 






9||/9' 



_n-q\\/q q\\/q 



nil 



(38) 



and case (xii) the resonances are equivalent to Eq. (plj). The eighth term of 'E.ijkh has a resonance of the type given in case (xi), 
and additional resonances at (xiii) ai6^„j — 6f„a2 = and (xiv) 036^^ — 6^j(Xi = 0. Repeating the procedure from above, we 
get for case (xiii) the solution 



hnqu/q 



± 



2m,f. 



hnqn/qj h?{n — qu/q) 



q\\/q 



^n . ^Ti 



(39) 



and in case (xiv) gives resonances equivalent to the result 
given in Eq. (B4). Again, when considering a two-level quan- 
tum well in our configuration, some restrictions apply to the 
quantum numbers. If we again compare Eqs. (Al) and (HO) 
we see that if pump field one (index k) is s-polarized, then 
V = n, and if it is p-polarized, then v ^ n. Additionally, if 
pump field two (index h) is s-polarized, m = I, and if it is p- 
polarized, ra j^l. This has the consequences that (i) the quan- 
tum numbers n and m can be chosen arbitrarily when both 
pump fields are s-polarized, (ii) when both pump fields are p- 
polarized we either get m ~ n and I ~ v, or we get m — v 
and I = n, (iii) when pump field one is s-polarized and the 
other one p-polarized we get either m — v 01 1 ~ v, and (iv) 
in the opposite case we get either m = v 01 m = n. These 
conditions are summarized in Table l|, and the contributions 
from Eqs. (p4|)-(P9|) to the resonances in Fig. y\ are shown in 
Table |l| for the valid combinations of quantum numbers. It 
should be noted that in Eq. (p4|), the combinations of quantum 



numbers that give rise to the resonances "b", "e", "i", "h", and 
"1" are going into resonance "m" after they have reached the 
line at iojioix — \. None of the other equations contributes to 
resonance "m". 



For the ninth term of ^ijku the resonances are at (xv) 
a2&Li - &L"2 = 0, (xvi) 046^^ - 6^,„a2 = 0, and (xvii) 
'^4^?m ~ bt„,a7 = 0. Af ter insertion of the relevant a's and 
6's from Eqs. (|Cl|)-( C13 ), fc|| = tiuj/cq and a normalization 



of (Jii to the vacuum wavenumber, we resulting second order 
equations can be solved with respect to w as a function of 
q\\/q. Then case (xv) is equivalent to Eq. (|33|), and cases (xvi) 
and (xvii) to Eq. (P5|). Finally, in the tenth term of S^fe/j a 
resonance of the type given by case (xvi) occur Two other 
resonances are located at (xviii) oifofj,,, + fo^/Oa = and at 



(xix) a^bt, + - 



,ai 



0, respectively. Inserting the a's and 



6's given in Appendix^ and using the same substitution and 
normalization as above, case (xviii) gives 



TABLE III. Resonances generated by Eqs. (p0|)-(|40|) are shown as a function of the valid combinations of quantum numbers (n, m, v, I) 
and the sign appearing in front of the square roots. In each of upper and lower parts of the table, the upper row shows the generating equation 
and the next four rows show the values of the quantum numbers, which can take the value 1 or 2 in a two-level quantum well. The last four 
rows show the resonances resulting from use of the quantum numbers in the respective equations for each sign + and — , the first two of these 
rows being associated with the normal term, and the last two with the interchanged term. A zero in the last four rows refers to oj = 0, and the 
letters a~l refers to the resonances shown in Fig. pi A "*" is used when the value of a quantum number is indifferent, and a "-" in the output 
field appears when the result is outside the shown range in Fig. u\ Since Eg. (t33p is a linear solution in cj the sign does not apply, and the result 
is listed under the sign "+" for simplicity. It should be noted that in Eq. (p4), that the combinations of quantum numbers that give rise to the 
resonances "b", "e", "i", "h", and "1" are going into resonance "m" after they have reached the line at uj/uj2i ~ 1- 
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niecl n + q\\/q 



± 



''e'-'O 



q\\/q 



hnq\\ /qn-q\\/q \ \ hnq\\^ /q n - q\\ jq 

and case (xix) has a solution equivalent to the one given in 
Eg. (pq). As before we find by a comparison of Eqs. ( |A11 ) and 
( pX]| ) that some selection rules appear when choosing a two- 
level quantum well in our configuration, since when pump 
field one (indexed fc) is s-polarized we get v = n, and when 
it is p-polarized, v ^ n. Similarly, when pump field two 
(indexed h) is s-polarized we get / = v, and when it is p- 
polarized, I ^ v. Then, if both pump fields are s-polarized we 
may in a two-level quantum well choose m — n or m ^ n. 
In the case where both pump fields are p-polarized, the re- 
sult is identically zero. In the case where pump field one is 
s-polarized and pump field two is p-polarized we may choose 
either m — I or m — v, while in the opposite case we may 
choose either m — I or m ^ n. As before, these condi- 
tions are summarized in Table O, and the contributions from 
Eq. (40) to the resonances in Fig. are shown in Table |l| for 
the valid combinations of quantum numbers. 

In the linear conductivity tensor [Eq. (^1 resonances occur 
when 026; 



2mecg 



4 

nm 



'^^bnln — 0' where 



1 ^Qn 



(41) 



The solutions are q\\ = or w = 0, independent of the val- 
ues of n and m. Adding this resonance to the ones we found 
in Eqs. (|30|)-(p^ the resonances appearing in Figs. ||-Hhave 
been identified. Q. E. D. 

While most of the resonances described above and shown 
in Fig. M are clearly pronounced in Figs, pj-pl the resonance 
named "m" does not appear so clearly, although in Figs, g 
and H the curves indicate that something is present around 
the position of "m". This resonance is striking by the fact 
that it approaches the Fermi wavenumber when the frequency 
approaches zero. It might also be appropriate here to men- 
tion that the resonances named "a" and "b" have the asymp- 
totic value of (7[| = 1/n in the low end of the normalized 
g||-a;-spectrum, and that the resonances named "c" and "d" 
approaches g|| — n for high values of gy /q and low values of 
(jj /u)2i- The resonance named "h" is the interband resonance. 



h?i 



m/q) 



q\\/q 



(40) 



understanding of the importance of the different resonances. 
However, since it would be rather difficult in an experiment 
to follow such a path, and since the exact positions of the 
resonances probably will be shifted in a practical situation, 
we have chosen not to do so. We have in stead in Figs. || 
^ and K^ plotted the intensity of the phase conjugated field 
along linear cuts in the normalized a;-g||-plane at q\\/q = 0.4, 
q\\/q — 3.0, and lo/lj2i = 1.5, respectively. Following the 
curves in Figs, pj-p^ along their respective path on Figs. H- 
g, the appearance and dissapearance of each resonance along 
the path is easily identified. From Figs. p[-p^ it also appears 
that some of the regions in Figs. ^ and Q with high density of 
isophotes are zeros rather than resonances. 



10'- 




UJ/UJ21 



V. DISCUSSION 

To give an impression of the magnitude of the phase conju- 
gated response, we have in Figs. |]-g| highlighted the isophotes 
with magnitude of 10^°m''AV^ and lO'^^m^/W^ by drawing 
them with a long-dashed curve and a short-dashed curve, 
respectively. Their positions in the normalized g||-a;-plane 
shows quite clearly that most of the area reachable within a 
single-mode experiment should produce a phase conjugated 
response of a magnitude comparable to what one gets from 
second-harmonic generation (compare Refs. pT|- p^p^|Jl 9| , p5| ) . 

Knowing the positions of the resonances in the normalized 
g||-a;-plane, one could of course be tempted to plot the magni- 
tude of the phase conjugated response along paths following 
each of the resonances (e.g., following the path of resonance 
"i", and its continuation into "m") in order to give an improved 



FIG. 8. The phase conjugation reflection coefficient is shown for 
the four combinations of polarization presented in Figs. §-H in the 
normalized angular frequency range 0.1 < ujIujix < 10 for a con- 
stant value of the parallel wavevector, qy = OAq. Thus the four 
curves represents the ppp (dash-dot curve), pps (fully drawn curve), 
ssp (dashed curve), and spp (dotted curve) configurations. 



One of the resonances are of special interest, namely the 
resonance at the interband transition frequency, which experi- 
mentally is rather easy to tune into. Until now, resoaant four- 
wave mixing has been studied in other contexts,ETl£3 but al- 
ways at the point {q^^ , oj) ~ (0, UJ21) in the q^^-us-plane. To go 
beyond that, we have plotted the phase conjugated response in 
the case where the interband transition is resonant (along the 
linear path in the normalized g||-a;-plane where oj = u!2i) in 
Fig-O 
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FIG. 9. The phase conjugation reflection coefficient is shown for the 
four combinations of polarization presented in Figs. B-lq in the nor- 
malized angular frequency range 0.1 < lo/lu2\ < 10 for a constant 
value of the parallel wavevector, gy = 3.0g. The ppp configuration 
result is drawn using a dash-dot type of curve, while the pps, ssp, and 
spp configurations are drawn using dotted, dashed, and fully drawn 
curves, respectively. 

In configurations with only a single source field in the field- 
matter interaction, such as, e.g., in linear response, second- 
harmonic generation, photon drag, and photoemission the so- 
called self-field approximation has proven to be quite effec- 
tive. The founding argument to use the self-field approxima- 
tion is that the dynamics across the quantum well (in the z- 
direction here) are dominating over motion in the plane of the 
quantum well (x-y-plane here). Let us as a test in the follow- 
ing look at the consequences of applying the self-field approx- 
imation in the present case of degenerate four-wave mixing, 
where three incident fields are present. 

Working within the framework of the self-field approxima- 
tion, we observe from Eq. ^ that the phase conjugated re- 
sponse would have been limited to the cases where nonlin- 
ear and linear current densities is produced in the z-direction. 
Hence, only tensor elements with i ^ z would contribute. 
Then, from Table || we observe that the contributions from (i) 
the two cases where the pump fields have the same polariza- 
tion and the probe field is s-polarized (sss and pps), and (ii) the 
mixed-pump configurations spp and psp would have been ne- 
glected. Thus, the data presented in Figs, g and ^ would have 
been absent. While this is certainly a good aproximation in 
the pure s-polarized case, the argument is not so good in cases 
with pump or probe dynamics in the z-direction. Using the 
argument of the dominating z-dynamics, it is striking that the 
mixed-pump configurations with s-polarized probe field sur- 
vives the self-field approximation while the two others do not, 
because we would expect more dynamices in the z-direction 
from the latter two. Another interesting conclusion is that with 
the loss of Fig. || we would also lose the resonances named "j", 
"k", and "1" in Fig. M. At the same time we would keep the 



essentially nonresonant ssp case. Comparing the raw arripli- 
tudes of the different configurations we can see from Figs. ^-^ 
and p|-|lO| that in most regions of the ^n -a;-plane, the ppp con- 
figuration gives a response that is a few orders of magnitude 
larger than the other configurations, but we also observe that 
the three other cases have resonances around q\\/q = 0, while 
the ppp-case do not. Thus, for near-normal incidense of the 
probe, the phase conjugated reflection coefficient is larger for 
some of the mixed modes than for the pure p-polarized con- 
figuration, indeed leaving room for experiments that cannot 
be described within the framework of the self-field approxi- 
mation. 



u 




0.1 0.1 



FIG. 10. The phase conjugation reflection coefficient is drawn on 
a logarithmic scale for the four combinations of polarization pre- 
sented in Figs. H-n in the normalized parallel wavevector range 
< Qw/q — kp/q for a constant value of the angular frequency, 
uj = 1.5a;2i. In the strip to the left, the abscissa is linear, while it is 
logarithmic in the right part of the figure. The scale of the ordinate 
is the same in both frames. The upper curve (dash-dot) shows the 
result for the ppp configuration of polarizations, while the dashed 
curve shows the pps result, the fully drawn curve shows the ssp re- 
sult, and the dotted curve shows the spp result. 

All in all, we may conclude from the above discussion that 
although the argument behind the self-field approximation re- 
mains intact, when one allows more than one incident field 
to participate in the interaction (as in, e.g., sum- and differ- 
ence frequency generation, or degenerate four-wave mixing), 
one should be careful in applying the self-field approximation 
in cases where mixed polarizations of the incident fields are 
allowed. 

Outside the resonances the influence of the relaxation time 
is insignificant, but around the resonances the choice of re- 
laxation time has a great influence on the width (in the q^\- 
space) and amplitude of each resonance. Choosing adequate 
relaxation times T„m is a difficult problem and it appears from 
Fig. O how big impact the relaxation time has on the phase 
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conjugation reflection coefficient. The intraband relaxation 
time in the occupied state (rn) has been chosen in accordance 
with Ref. ^to be 3fs. For the unoccupied state the relaxation 
time T22 (see Fig. g) has been chosen to approach infinity. In 
the present case where also interband transitions contribute to 
the phase conjugated response, the intraband relaxation time 
is of little importance, and thus it is the choice of interband 
relaxation times (here T21 and T12) that are critical. In the 
present calculation we assume no relaxation from state |1) to 
state 1 2), letting T12 — > 00. 



Is 

o 



10'^ - 



10" 



10"= - 



I 
/ 

'V') 
'■■' ' 
r\ 

_/ 



0.1 




1000 



q\\/q kp/q 

FIG. 11. The phase conjugation reflection coefficient is shown on 
a logarithmic scale for the four combinations of polarization pre- 
sented in Figs, g-g in the normalized parallel wavevector range 
< 911/'? — kp/q when the value of the angular frequency is ex- 
actly equal to the interband transition frequency, ui — uj2i- As in 
Fig. hd the strip to the left shows the range < q\\/q < 0.1 with 
linear abscissa, while the rest is plotted with logarithmic abscissa. 
The scale of the ordinate is the same for both frames. In this figure, 
the dash-dot curve corresponds to the ppp case as in the previous 
figure, but the dotted curve to the pps result. The fully drawn curve 
corresponds to the ssp case as before, and the dashed curve to the 
spp result. 

The phase conjugation reflection coefficent has in Fig. |lj 
been plotted for four values of the relaxation time from state 
|2) to state |1), namely (i) 30fs and (ii) 20pfe, which are typ- 
ical values one would find for bulk copperEil at (i) room tem- 
perature and (ii) at 77K, (iii) 3fs, and (iv) 2ps. The value in 
case (iii) is obtained by a conjecture based on the|^fference 
between measured data for a lead quantum wellE£l and the 
bulk value for lead at room temperature. The difference be- 
tween the_celaxation times measured by Jalochowski, Strozak, 
and ZdybL^ is for two monolayers approximately one order of 
magnitude. Case (iv) is included to see the effect of raising the 
value of the relaxation time one order of magnitude, thus es- 
sentially assuming a better conductance than in case (ii). The 
values (i)-(ui) are the same values as we chose in our descrip- 
tion of the single-level quantum- well case where only intra- 



band transitions were allowed£3 but since the interband tran- 
sition is of a more bulk-like character we have in the present 
calculations chosen T21 — 200fs. We notice that in the case 
where both pump fields are s-polarized (polarized in the plane 
of the quantum well), the phase conjugated response does not 
vary as a function of the interband relaxation time, whereas 
in the other three cases the general tendency is that they have 
larger magnitudes for larger values of the relaxation time. 




10-24 H 



FIG. 12. The phase conjugation reflection coefficient is shown for 
interband transition resonance for different values of the interband 
relaxation time T21 € {3, 30, 200} femtoseconds, and 2 picosec- 
onds. The fully drawn curve corresponds to 200fs, the long-dashed 
curve to 30fs, the short-dashed curve to 3fs, and the dotted curve to 
2ps. 



VI. CONCLUSIONS 

Our main conclusion from this work is that DFWM in a thin 
metallic film gives rise to several resonance structures even in 
the propagating regime of the (/n -spectrum. Furthermore the 
coupling by the phase conjugation reflection coefficient is of a 
magnitude that is well within experimental reach. Thus, also 
single mode excitation in the experimentally feasible regime 
(up to arowid n = 3) should-^e possible by use of the stan- 
dard OttoO or Kretschmann23 techniques, and a qualitative 
comparison with the present work should be possible. How- 
ever, for a better quantitative comparison in a specific system, 
it will be necessary to refine the numerical calculation by, e.g., 
abandoning the IB model in favor of one of the flavors of the 
KKR, LAPW or LMTO models, although such a task may 
prove to be strenuous. 
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APPENDIX A: NONLINEAR CONDUCTIVITY TENSOR 



Under the assumption that the electron dynamics is free-electron-like in the plane of the quantum well the nonlinear response 
function S(z, z'; (jj| , fcy , uj) [given by Eq. (m, and with tensor elements '^tjkh] can be obtained from the results established for 
E!'^(z, z', z", z"'\ q\\ , fc|| , Lu) in Ref. ^ Upon integration over z" and z'" [and use of Eq. (g)] one gets 

^ijkhiz,z';q\\,k\i,uj) 



8^3 (27r)2 (iuj)3 ^ 



i0„rn{K\\ +q\\,K\\) -to 



fn{K\\ + q\\) - fv{K\\)\ 1 



f fiJH -fc||) -/f(^||) 

OJvm (k|| , K|| ) - 2w y (i^; (k|| , K jl - fc|| ) - W 



^nv (K|| + Qll , Kll ) + tJ y Cjr,l (K|| + (TJI , K|| - fc|| ) _ 

X / jh,mi{z'";2K\\ -k\\)dz"' I jfe,;„(z"; 2k|| - fc||)dz"jj\«„(2;'; 2k|| + q\\) 
fi{K\\ - k\\) - f„i{Ki\) //(k|| - A:||) -/„(kj| -fc|| +g||)\ 



+ 



y(I;im(K|| - fcl[,K||) - W t;>i,((K|| - fc|| +q\\,K\\ - fc||) + wy Ulvrn(K\\ - fc|| + ^jl , '^H ) 

/((^ll - fc||) -/^'(^ll -^11 +'?i|) , fni^w + q\\) ~ fv{K\i ~ kii + q\\)\ 1 



ycJu;(Ki| - fc|| +<7||,K|| - fc||) + UJ a)„„(K|| +(?j|,K|| - fc|| +q\\) - uj ) cj„;(k|| +gj|,Kji - A:||) 
X / 3h.mi{z"' \2k\\ ~k\\)dz"' j jfe,„„(z";2K|| ~ k\\ + 2q\\)dz" jj^i^{z' \2k\\ - 2k\\ + q\\) 

f M^W +gi|)-/m(At||) /;(^|| + g||) - fv{K\\ - fcy + g||) \ 1 

y Wim(K|| +g||,K||) +0; a)„;(K|| - fc[| + gji , K[| + (7J| ) - t^; y ^^^^(kii - fc|| +gj|,K||) 
/K^ll +9i|) -/u('2|| -^11 +'?i|) , ./n(«|| +9i|) -/«('«|| -^11 +9i|) \ 1 



^(I)i,;(kII -fc|| + (7|| , K|| + 911 ) - tj iJj„v{k\i +gJ|,K|| -fc|| +g||)-cjy i^nilKll +9ih^|| +9||) -2w^ 
X / jh,iv{z"' ■■,2k\\ - fc|| +2gj|)dz"' / jfc^„„(z"; 2k|| - fc|| + 2(7||)dz"jj,,„/(z'; 2k|| +g||) I ji,„„(z; 2k|| +q||)d^K||. (Al) 



APPENDIX B: ON THE SOLUTION TO THE INTEGRALS 
OVER K I IN THE LOW-TEMPERATURE LIMIT 

In this appendix we discuss how analytical solutions to 
the integrals over the electronic wavevector, Ky, appearing in 
the linear and nonlinear conductivity tensor may be obtained, 
and for simplicity the discussion is limited to cover the low- 
temperature limit. These integrals can, when scattering takes 
place in the x-z-plane, be expressed as a sum over terms of 
the general type 



T^Jn,{a],{b},s) = 



OO /"OO 



OO ■/ — CXD 



Uk^li'^kKx + bk] 



number n, which is a positive nonzero integer, (ii) a set of 
real quantities, {a} = {ai, . . . , a^} appearing in front of the 
integration variable k^ in the denominator, (iii) a set of com- 
plex nonzero quantities, {b} = {bi, . . . , 6^} appearing also in 
the denominator, and (iv) the real quantity s representing the 
displacement (in the a;-direction) of the center of the Fermi- 
Dirac distribution function from {k^, Ky) = (0, 0). The quan- 
tity s together with each element in the set {a} are in general 
functions of the parallel components of the probe and pump 
wavevectors, q^\ and fc||. Each element in the set {b} is fur- 
thermore a function of r, the relaxation time. 



dfXi^ (Jjf%1j . 



(Bl) 



where p,k, (3 are nonnegative integers, and q is an even non- 
negative integer The functions depends on (i) the quantum 



The combinations of p and q needed in Eq. (Bl ) in order to 
solve the integrals over Ky in the nonlinear conductivity tensor 
are {p,q) G {(0,0), (0,2), (0,4), (1,0), (1,2), (2,0), (2,2), 
(3,0), (4,0)}, and /3 e {1,2,3}. However, functions with 
(3 = 2 and /3 = 3 can be expressed in terms of functions with 
/3 = 1 in the following way: 



,Fp(n, 01,02,61,62,5) = 



aiTl Jn, ai, 5i, s) - a2Tl {n, 02, &2, s) 



:r'^ (71,01,02,03,61,62,63,5) 



0162 - 0261 
alJ^lq{n,ai,bi,s) 



(B2) 



alJ^L{n,a2,b2,s) 



+ 



0'lHq{'n,a3,b3,s) 



(0261 - 620i)(a36i - 6301) (0261 - 620i)(a362 - 6302) (0361 - 63ai)(o362 - 6302) 

(B3) 
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Eqs. ( [B2| ) and ( |B3[ ) are given with the provision that the val- 
ues of the different Ofc are nonzero, k G {1, 2} in Eq. (B2) 
and k G {1,2,3} in Eq. (B3). If any afc, for instance 
oi, becomes zero, we see from Eq. (Bl) that the order 
(in Kx) of the denominator becomes smaller by one. This 
implies that :r2^(n,0,a2, 61,62,5) = Tl^(n,a2,b2,s)/hi 



in Eq. (82). The similar 



•^p<7('^, 0,02, as, 61, 62, 63,3) 



pq 

conclusion for 

-I pq 



Eq. ( |B3| ) is 

(n, 02,03, 62,63, s)/6i. 
Analogous reductions applies for any other Ofc = 0. 

In the low-temperature limit the Fermi-Dirac distribution 
function is zero outside the Fermi sphere and equal to one in- 
side, and it is therefore advantageous to shift n^ by — s, and af- 
terwards carry out the integrations in polar (r, 9) coordinates 
Using Kx — r cos 9, Ky — r sin 
integrals to be solved are of the type 



and duxd-Ky = rd9dr, the 



Tl {n,a,b,s) 



a{n) 
"'0 



r(r COS — s)P(r sin 0)' 



as 



ar cos ( 



d9dr, 
(B4) 



dropping the now superfluous index on a and 6. The upper 
limit of the radial integration is a(n) = \/kp — [nn/dY, 
kp > Tm/d. If kp < TTTi/d, the Fermi-Dirac distribution 
function is zero, and thus the integral vanishes. Physically, 
a{n) may be characterized as the two-dimensional Fermi 
wavenumber for electrons in subband n. 

Since the following treatment is a formal solution of 
Eq. (B4), we will also drop the reference to n for brevity, 
letting a = a{n). To solve Eq. (B4), let us make the sub- 
stitutions 



V 



aa 



au, 



(B5) 



and thereby turn Eq. (B4) into 

^pq{a,V,s) = 



JO 



^"^ u«+i(aucos6' - s)P{l- cos2 



^9/2 



T] + u COS 9 



-d9du, 
(B6) 



i.e., compared to the possible values of p and q, an expression 
where the angular integral is expressed as a sum of terms of 
the form cos'' 9 in the nominator, where /i G {0, 1, 2, 3, 4}. To 
carry out the angular integrals we put t = exp{i9) so that the 
integrals become of the type 



2ir 



cos''6' 



-J9 = 



1 



V 



ucoat 



2\h 



(1+t^) 



2''to/ t'^{t-t+){t-t^) 



dt. (B7) 



In Eq. (B7), the poles at t± in the f -plane are located at 



t± 



± 



(B8) 



and the integration runs along the unit circle. Since <+t_ = 1, 
one of these poles is located inside the unit circle while the 
other is outside. When h > Q there is an additional pole of 
order /i at t = 0. Using the unit circle as contour, residue 
calculations give the nontrivial solutions 



27r 



d9 



2ti 



ry + ucos^ ^rf - v? ' 



COS0 ,^ 27r 

-rft^ = — 

77 + u cos 9 u 



^/rf' — u^ 



cos^ 9 ,^ 27r77 

-d9 — 7T- 



rj + u cos ( 



■\/ry^ — u^ 



d9 = - -\ ^ 



T] + u cos ( 

^" cos^ 9 
r] + u cos ( 



u U"^ 



\nf 



■d9^ 



2T:Tf 



vV 



- 1 



(B9) 
(BIO) 
(BID 
(B12) 



- ^. (B13) 



To finish the formal solution, (i) insert these results into 
Eq. (B6), (ii) carry out the elementary radial integrations (see, 
e.g., Ref. p^. Sec. 2.27), (iii) backsubstitute rj, and (iv) check 
convergence for a ^ 0. Step (iv) can be carried out by use 
of a binomial series expansion of the square roots appearing, 
and a comparison the the result one gets by setting o = al- 
ready in Eq. (B4). The solution to the integrals appearing in 
Eqs. (Q) and (kl) are then found in a straightforward man- 
ner, but since the algebraic expressions are rather long, we 
will omit presenting them here [they can be found in Ref. p^ 
together with explicit expressions for the case where o = 0]. 



APPENDIX C: DENOMINATOR COEFFICIENTS 



Ol 



0.2 



hk\\ 

TOe 

me 

h 



as = (<?|| +fcj|), 



04 



(9[| -fc||), 






6^ = 



6^ = 



6* = 



6^ = 



n 



h 



h 



h 



hkj 

2me 
2me 

H_ 

2me 

H 

2me 

h 



-UJ -1T„ 



-U! - 1T„ 



-UJ -1T„ 



^(e« -em) + :7^(g|| -/e||) " ^^' 



6« 






6« = 



h 
1 

h 
1 



2me 

h 

2mc 

!m 

2me 
2m^ 

27Tl« 



2 ;_-l 
nm 7 



(CI) 

(C2) 
(C3) 
(C4) 

(C5) 

(C6) 

(C7) 

(C8) 
(C9) 
(CIO) 



.(g„-2fc||)+^-iT-^, (Cll) 
.(29||-/ci|)-c^-iT„"„\, (C12) 

{kn-2qn)-i,-iT-^. (C13) 
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